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Abstract 

The lifetime of a stored beam due to the Touschek effect is calculated 
for arbitrary ratios of beam height to beam width. A variation of 
the beam envelopes is taken into account, i. e. the derivatives of 
the horizontal and vertical amplitude functions and dispersions are 
included. The calculation is done for arbitrary energies in the rest 
frame of the colliding particles. 



1. Introduction 

Coulomb scattering of charged particles in a stored beam causes an exchange 
of energies between the transverse and longitudinal motion. It changes, therefore, 
the betatron and synchrotron oscillation coordinates of the colliding particles. One 
consequence is the Touschek effect which is the transformation of a small transverse 
momentum into a large longitudinal momentum due to the scattering. Then both 
scattered particles are lost, one with too much and one with too little energy. The 
amplification of the momentum change is a relativistic effect so that the change of 
the longitudinal momentum is increased by the Lorentz factor 7. 

The Touschek effect is different from intra-beam scattering which is also caused 
by Coulomb scattering. The intra-beam scattering, however, is a multiple scattering 
which leads to diffusion in all three directions and, primarily, changes the beam 
dimensions. The Touschek effect, on the other hand, is a single scattering effect 
which leads to the immediate loss of the colliding particles. Here only the energy 
transfer from the transverse to the longitudinal direction plays a role 

The Touschek effect was observed for the first time in the small electron storage 
ring ADA [1] and it turned out later that it is a serious problem in storage rings 
at low energies and in synchrotron light sources where it can reduce the lifetime 
considerably. Theoretical investigations were made in [1, 2] for the non-relativistic 
case, where the velocities of two colliding particles are non-relativistic in their center of 
mass system. The theory was then extended in [3] to the ultra-relativistic case, and in 
[4] to arbitrary energies. Dispersion was included in [5]. In all these calculations only 
the transfer of horizontal momentum into longitudinal momentum was considered 
whereas the transfer of the vertical momentum was neglected. In [6] a 100% coupling 
of horizontal and vertical betatron oscillations, i. e. the special case of a round beam, 
was considered for the non-relativistic case. 

The main aim of the following calculation is to take into account both the hori- 
zontal and the vertical betatron oscillations. In that case one has a two-dimensional 
distribution of the transverse momenta which leads to a weaker dependence on the 
maximum stable energy deviation than a one-dimensional distribution. This can 
increase the Touschek lifetime by a factor of two as compared to the flat beam ap- 
proximation at energies around 1 GeV, but the correction can be larger at lower 
energies. 

At the same time a variation of the beam envelopes is taken into account. The 
derivatives of the amplitude functions (5 X and (5 Z and of the dispersions D x and D z 
increase the angle between colliding particles. But the derivatives of the amplitude 
functions play a role only at positions where the dispersion does not vanish, or to 
be exact, where f3 x>z D' x z — P' xz D XjZ /2 is different from zero. In that case they can 
increase the loss rate and decrease the Touschek lifetime. 

We will use here the same method as for the calculation of the intra-beam scat- 
tering [7] . Thus we will consider the collision of two arbitrary particles in their center 
of mass system using the complete M0ller scattering cross section, i. e. assuming ar- 
bitrary energies in the center of mass system. After transforming into the laboratory 



2 



system we will calculate the number of collisions with those momentum changes which 
lead to the loss of both particles. When averaging over the positions and angles of 
the colliding particles we will assume Gaussian distributions for all coordinates. 



2. Change of the momenta of two colliding particles 

The momenta p\ and p 2 of the two colliding particles before the collision are given 
in the laboratory coordinate system {s,x, z} by the two vectors 




Pl,2 = Pxl,2 \ (1) 



where s indicates the longitudinal, x the horizontal, and z the vertical direction. We 
define a coordinate system with the axes j, k, and / which are parallel to p\ + p 2 , 
Pi x p 2 , and (pi + p 2 ) x (pi x p 2 ), respectively. The two momenta then take the form 




Pl,2 = Pl,2 | | (2) 

j,k,l 

where Xi an d X2 are the angles between the vector pi +p 2 and the vectors pi and p 2 , 
respectively. They are given by the two relations 

P1P2 = PiP2Cos(xi + X2) (3) 

and 

Pisinxi = P2SUIX2 (4) 

When we apply a Lorentz transformation parallel to p\ + p 2 we obtain for the 
momenta p x and p 2 in the c. o. m. system the representation 

/ 7t(cosxi,2 - Pt/M \ 

Pl,2 = Pl,2 

V ±sinxi,2 J m 

The transformed energies E\ and E 2 are given by (see App. Al) 

£1,2 = 7^1,2(1- PtPi,2 cos X i >2 ) = (E 1 +E 2 )/(2^ t ) (6) 

where f3 t is the relative velocity of the c. o. m. system, ^ t is the Lorentz factor of the 
transformation, (3 lj2 are the relative velocities of the two particles in the laboratory 
system, and the bars denote all quantities in the c. o. m. system. (3 t is determined by 
the condition that the sum of the two momenta vanishes in the c. o. m. system and 
is given by: 

a = + p 2 | = Pi 7i cosxi + P2I2 COSX2 (7) 
E! + E 2 " ' " ' 71 + 72 
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We may now assume that the quantities {jp\— P2) 2 /p 2 ^ v\\,il V\,ii an d Pli^/Pi^ 
and therefore also xl,2i are small as compared to one, and Eqs.(3) and (4) become, if 
we replace p sl:2 by Pi,2-(p 2 i, 2 +P 2 i,2)/( 2 Pi,2): 



(Pxl/Pl-Px2/P2) 2 + (Pzl/Pl-Pz2/P2) 2 = (Xl + X2Y 



and 



P1X1 = P2X2 
Finally one obtains with xi ~ X2 ~ X 

4X 2 = (Pzl -Px2f/P 2 + (Pzl -Pz2) 2 /P 2 



(8) 



where p is the mean value of all momenta in the bunch. Eqs.(5) and (6) then simplify 
to (see App. Al) 

^1 + 7 2 X 2 /2 \ 




Pl,2 = ±P 



/ 



j,k,l 



and 
with 



E 1;2 = E/~f t 



(9) 



(10) 



Pi -P2 

IP 



e = 



Pxl - Px2 
P 



C = 



Pzl - Pz2 
P 



4 X Z = e 2 + ( 2 



We now assume that the momenta are almost perpendicular to the j-axis in the 
c. o. m. system, which means (Eq.(9)) 

e 2 (l + 7 Y) « Y or e « (4- 7 2 e 2 )x 2 ~ 4 X 2 
which can be written as 

(P1-P2) 2 < 7 2 ((Pxl-Px2) 2 + (Pzl-P 2 2) 2 ) 

This assumption gives 



and one obtains 



\Pi\ ~ IP2I ~P 
/ 



Pl,2 = P 







V ± ^ / j,k,i 

The relative velocity (3 t of the c. o. m. system simplifies to 

Pt = P cos x 

where (3 is given by the mean value p of all momenta and j t is 



It 



r 



1 _ 02 ! + ^2 7 2 X 2 
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11) 



(12) 



(13) 



(14) 



After the collision the absolute values of the two momenta are not changed in the 
c. o. m. system and the new directions of the momenta can be described with help of 
the two angles -0 and 0, where -0 is the angle between the momentum of the scattered 
particle p x and the (longitudinal) j-axis, and is the angle between the projection of 
p-y on the £-/c-plane and the /c-axis. The momenta after the collision are then given 
by 



Pl,2 



±P 



( COS 

sin -0 cos 
\ sin ip sin 



±p sin x 



3,k,l 



( cos-0 

sin -0 cos 
\ sin if) sin 



(15) 



3,k,l 



The inverse Lorentz transformation gives the two rotated momenta in the laboratory 
system: 



Pl',2 



lt {±p' +p t E/c) 

±P' k 
±P'i 



3,k,l 



±7t sin x cos -0 + cos x 
p | ± sin x sin -0 cos 
± sin x sin -0 sin 



The change of the momenta due to the collision is 



(16) 



±p sin x 



I 7 t cos -0 

sin-0 cos0 
y sin -0 sin — 1 



±PX 



3,k,l 



I 7 t cos -0 

sin-0 cos0 
y sin -0 sin 0—1 



3,k,i 



(17) 



With = Jpf 12 + pf 12 + pg 2 « 2 + (pjg j2 + p/ 1 2 2 )/(2pj 12 ) we get finally 



Pl,2 - Pl,2 ~ - Pjl,2 



± p 7 t x cos -0 



(18) 



3. Scattering cross section 

The probability for scattering of one of the two colliding particles into the solid 
angle dCl is given by the M0ller scattering cross-section in the c. o. m. system [8]: 

d ° = t4 (i 1 + i) 2 (-^ - -4s) + -4^ + d« (19) 

4^2 ^2/ V sm 4$ snr<3v sin 2 $ / 

r p is the classical particle radius, 7 is 7/7*, and $ is the angle between the momenta 

— 

before and after the collision, i. e. between the £-axis and p 1 . In order to determine 
all collisions with an energy change larger than the maximum stable energy deviation 
we have to integrate over 05 and with the conditions 

O<0<2vr, < < m 

where -0 m is given by the maximum stable momentum deviation Ap m . If one of the 
two particles is scattered into the region n — -0 m < < n the other one is scattered 
into < < m so that both regions are included. With Eq.(18) we obtain for Ap m 
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Ap m = \p{ t2 - Pl,2\m = PX.1t COS l\) Tl 

or 



T 5 m Wj + /3 2 7 V . , fon , 

COSIpm = = < 1 (20) 

ItX IX 



with 

Eq.(20) gives also the condition 



Apr, 
P 



r) 2 f) 2 

2 > 2 = <>m V (21) 

A — Am _,2^1 N -2 V / 

7 2 (1-/^J Y 

With <if2 = sin tp d<j) dvjj one gets for the total cross section 

r 2 ri>m /-27T / . 1 \ 2 / 4 3 \ 4 \ _ _ _ 

^ = 74 / / + + + 1 sin^## (22) 

47 2 Jo 7o V v sm 4 $ siirSv snT$ / 

With the relation 

p' e = p cos $ = p sin ijj cos 
which follows from the definitions of the angles $, ip, and 0, one obtains 

sin 2 5> = sin 2 + cos 2 -0 cos 2 

Substituting tan0 = | cos-0| x tan it and integrating with respect to u and ip gives 

nr l ^ m ((, 1 \2/2(l + cos 2 ^) 3 \ 4 \ 

a = rf/ + (- 3T^ r) + T + l smutty 

27 JO \ v p z/ x cos 6 ip cosip / cos ip J 

2 7 2 

where (3 is given by 



? = ^ = 7(W = M ^1 (24) 

1 + p 7 x 1 + 7 X 

since x 2 ^ 1- The cross section a is parallel to p\ + P2, i- e. parallel to the £-axis, 
and it is transformed into the laboratory system by: 

a = - (25) 
It 
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4. Loss rate 



We want to calculate the number of collisions per unit time which lead to the loss 
of both particles. We shall integrate, therefore, the cross section with respect to all 
positions and angles of the two colliding particles satisfying the condition Eq.(21). 
The distribution of the positions and angles of electrons and positrons is, in good 
approximation, a Gaussian distribution: 



D , , , . 13 X I3 Z t x 2 + (a x x p + p x xf p ) 2 zj + (a z Z[3 + (3 z z' p ) 2 , 
Pp{x p ,x p ,z p ,z p ) - ex P{ } 

(26) 

with a x>z = —j3' XtZ /2. a x p and a z p are the standard deviations for the horizontal and 
vertical betatron distribution. The distribution of the synchrotron coordinates As 
and Ap is given by 



where o v and o s determine the relative momentum spread and the bunch length, 
respectively. 

The number of scattering events per unit time for a single particle moving with an 
angle of 1\ with respect to the momentum of the opposing particle is (3 re icP Q a where 
the relative velocity j3 re ic between the two colliding particles is 2/3csinx ps 2/3c%. P 
is the spatial density in the laboratory system. The total number of scattering events 
per unit time is obtained by integrating over all positions and angles of all particles 
with the condition Eq.(21): 

R = 2[3c J Pa X dV (28) 

P and dV are given by 

P = N 2 P s (A Sl , A Pl ) P s (As 2 , Ap 2 ) Pp(x pi , x' pi , z pi , z' pi ) P p (x p2 , x' p2 , z p2 , z' p2 ) (29) 
and 

dV = dAsi dxpi dzpi dApi dAp 2 dx' pi dx'^ 2 dz' pi dz'^ 2 (30) 

N p is the number of particles per bunch. When averaging the position coordinates 
over the whole beam we assume that the two colliding particles have always the same 
position, i. e. we take into account the following conditions: 

A A n A ^ n A P 2 n A ^ n A P 2 

As 1 = As 2 , Xfi X + D x =xp2 + D x , zpi + D Z = Zf32 + D z 



p p p p 

The integration in Eq.(28) is carried out in Appendix A2 and yields (r = (3 2/ y 2 x 2 ) 
r 2 cP x (3 z a h N 2 l^ (T /T m s 
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I is the modified Bessel function and the other quantities are given by 

1 1 Dl + Dl Dl + Dl 



2 — — (°l a % + °Wx(3 ~ °lf3° 2 z p) ( 32 ) 



Bl - w^r^ 1 - ^ + w^r^ - (33) 



B 2 , = - - — (1 - 

Bl- 



??35?l v '" w ') (34) 



4 n- 4 

r m = P 2 5 2 m (35) 
In order to simplify the representation we have introduced 

D x ,z = otx, z Dx,z + Px,zD' xz (36) 

and 

~< z = <z + oplz = °%, z p + °l{Dl, z + Dlz) (37) 



6. Touschek lifetime 

The number of particles lost per unit time is given by 

dJ ^ = -R = -aNl (38) 

The fact that always two particles are lost is taken into account by averaging over 
both particles, 1 and 2, in Eqs.(28) to (30), so that there are always two contributions 
to the scattering rate which differ only by the indices 1 and 2. Integration of Eq.(38) 
with respect to t gives 

1 

= —at + const. 

N p 

or 

1 N 

n p = = ; l° T (39) 

at — const. l+i\ at 

where N is the number of particles at the time t = 0. A lifetime can be defined 
by 

k - - (id < 40 > 
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after that the number of particles drops to half the initial number. The brackets 
denote the average over the whole circumference of the storage ring. Using the same 
convention as other calculations we write 

^ = ( ; T ' 2pCNp — F(r m , Bi, B 2 )\ (41) 

with 

F(r m , Bl ,B 2 ) = ^{Bl-BD r m /°°(( 2 + -f \^ - l) + 1 - 

-i(4 + i)ln^)e-U B ,)^ (4 2 ) 

where B\ — B\ is given by Eq.(34). A faster numerical integration is achieved by 
substituting r = tan 2 k , r m — tan 2 K m : 

F(r m ,B 1 ,B 2 ) = 2 V^CB 2 - B\) r m ((2r + 1) 2 - 1) /r + r - ^ (1 +r) 

- (2 + 1) In ^) e- B -/ D (5 2 r) rfK 

Eqs.(41) and (42) describe the most general case with respect to the horizontal and 
vertical betatron oscillation, the horizontal and vertical dispersion, and the derivatives 
of the amplitude functions and dispersions. Special cases with some simplifications 
will be discussed in the following sections. 

Figures 1 to 4 show F(r m , B±, B 2 ) as a function of r m = f3 2 5l l 5^ for different 
ratios of beam height to width and for different energies. 



7. Special cases 
7.1 Plane orbit 

In case of a plane orbit and without coupling of the horizontal dispersion one 
obtains 

D z = D z = (43) 

and &h is given by (Eq.(32)) 



1 



+ 



Dl + D 



'x/3 



With Eqs.(33) and (34) and with a z p = a z B\ and B 2 simplify to 

I R2„2 R2 



Bifi 



2/3 2 7 2 



n 2 ^2 n 2 



(44) 
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and the lifetime is given by 



^ = (^4^ F(r m , Bl ,B 2 )) (45) 



7.2 Flat beam 

For flat beams we assume 



^ » ^ (46) 
2 ^ pa ^ I 



PI PI 



and D z = D z = 0. The last condition (plane orbit) is not necessary but simplifies the 
calculation. Thus we obtain with Eq.(44) 



P, 



and 



B 2 a 2 

B 1 -B 2 - Px x 



and with Eq.(34) and ah = a p a x p/a x 

_ — PxPz a x 



P'l^l^h 



The argument of the Bessel function at the lower integration limit (Eq.(42)) is (x m = 

Sm/l) 

5 2 B 2 a 2 B 2 5 2 B 2 



2Y 



(T 'x(l (T 'x a 1f3 



2 7 V 2 



z(3 



If this value is large as compared to one we have 7 2 cr^g/ LB 2 <C 5^ . This means 
that the change •jtCzp/Pz of the longitudinal momentum due to the initial vertical 
momentum a z p/f3 z is smaller than the maximum stable momentum deviation 5 m so 
that the vertical momentum does not contribute to the loss rate. The Bessel function 
can then be written approximately as [11] 

e-^ I (tB 2 ) « —L=e-«*-V = ^Sf ex P {--^^} (47) 

The lifetime follows with Eqs.(31) and (47) to 

1 / r 2 c(B x N p r°°f f U 2 ( rjr m x 



To 



r((2+-Y(^-i)+i- vi+r 
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and F(r m , B u B 2 ) (Eq.(42)) is given by 



F(r m , Bl ,B 2) = .M# Tm rf( 2 + I) 2 (^- l) + 1 VTT 



2B, ./.,. V.v M • v 5 /r„, 



2r^ t> 1 + tJ y/TTr 
l\2/r/V m \ VI + t 



-(4+i)ln^exp{-^}^ (49) 



where e m is a generalized expression already used in other investigations: 

5 2 8 2 a 2 5 2 8 2 a 2 



(50) 



7.3 Round beam 

Here we assume again D z = D z = 0. Then we can define the round beam by 

VxP&xP Vz/3 

and Bi and B 2 are given by 

g _ Px a x _ £m g _ g 

F(r m , Bi, B 2 ) simplifies to 



(51) 



F(r m ,5 1? 0) = v^Bxr™ /°°((2 + i) 2 (^ - l) + 1 - 

A m V v vi + t * i T / T 



--(4 + -)ln^exp{-r5 1 }^L (52) 



With help of the integral [11] 



/ e~ pr . = exp{-ap/2}K D (ap/2) (53) 

Jo Jr(r + a) 



where K Q is a Bessel function, the integral in Eq.(52) can be solved exactly except 
for the terms with ln(r/(l + r)). Then one gets for small r m : 

F(r m ,B 1 ,0) = v ^B 1 (2 + B 1 )(K 1 -K o )e B ^ 2 + 0(^) (54) 
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where K 1 and K Q are Bessel functions with the argument B 1 /2. 

Thus for a round beam F becomes constant with decreasing r m . This is different in 
the case of a flat beam where F increases continuously with decreasing r m . The reason 
for the different behaviour is the different distribution of the transverse momentum. 
When the maximum stable energy deviation 5 m is decreased by — dS m more particles 
will be lost having the smaller transverse momentum difference Ap±—dAp±. In a one- 
dimensional distribution the number of particles which can be scattered additionally 
is proportional to 2dAp x , but in a two-dimensional distribution it is proportional to 
2irAp_ L dAp_ L with A 2 pj_ = A 2 p x + A 2 p z . Thus with decreasing Ap±_ the increase in 
loss rate is larger in a one-dimensional distribution, i. e. in a flat beam. 



7.4 Non-relativistic case 

The non-relativistic case is defined by f3 2 <C 1, i. e. (Eq.(24)) 

PVx 2 « 1 (55) 
which means also 7 2 x 2 <C 1 since x 2 <^ 1- 4>m is given by 

7 m &m 

cosf m = = — 

ItX IX 

since 7* = 7 (Eq.(14)). One obtains for the cross sections with Eq.(23) 

- " r i /Vx- 2 1 + ln M (56) 



and with x 2 = t/(P 2 7 2 ^ 



2^ 4 7 4 X 4 V SI IX, 



" = ^b^- 1 + ln 7?J (57) 



With Eqs.(40), (55) and (A2.5) one gets for the lifetime 

1 / cr 2 (3 x f3 z a h N p r°° / t 1 r \ _ rBl dr \ 



and 

F(r m , Sx, 5 2 ) = ^ (S? - S|) r m f° (— - 1 - \ In — ) e"^ 1 1 (rB 2 ) A (59) 

One can obtain this result also from Eq.(42) by neglecting r as compared one. Using 
the same notation as [9] one can write with D z = D z = 
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with 



Gio — 



Bio = 



2 2 2 
a z(3 a x _|_ a xf3 



7.5 Non-relativistic case for a flat beam 

For the non-relativistic case one obtains with Eqs.(47) and (58) and with j3 2 
7 t V = X = vWt), and D z = D z = 



J_ / r 2 p cfi x N p roo/r 1 \ n Tm \ c:i f re ™) dr ' 

Te \&nf37 3 (T x cr x pcr z p<7 s Jr m \r m ' 2 r) ^ r m >r 2 



and 



F(r m ,B 1 ,B 2 ) 



oo / r 1 T 

1 in — 

Tm \ 7~m 2 T m 



jexpj-— - ) 



Following existing representations (see f. e. [9, 10]) we can write (r = er m /e m , e m / 



i 

Te 



£f?' N * , Hi— - 1 - - In— ) exp { - e} "4 



87r/3 3 7 5 (7^t7 2/ 3(7 s ^e m 
87r/3 3 -f 3 a x(3 a z/3 a s a x Sl 



C(e m ) 
C(e m ) 



with 



C(e m ) 



--l--ln-)e-^ 



or 



— e 
2 



C(e r 



+ 



1 + ^i + ^ln— ) e" e - 
2 2 ej e 



filO zfi 



F(r m ,B 1 ,B 2 ) 



For D x = and a x = cr x this representation was obtained in [9]. 

7.6 Ultra-relativistic case 

The ultra-relativistic case is defined by 

(3 2 1 2 X 2 » 1 

which means (3 2 « 1 since x 2 <^ 1- With Eqs.(14) and (20) one obtains 



It = -, 

X 



cos ip. 
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The cross section follows from Eq.(23) with (3 2 = 1 — I/7 2 = 1 — 1/(7 x ) ~ 1 to 



2ttt* 



With Eqs.(25), (40) and (A2.5) and with (3 = 1 one obtains 



(65) 



!r = / ^f^fv r eM _ TBl}Io{rB2)dT \ (66) 



Using a definite integral [11] we can write 



^7" = ( o 7= 4 2 — 2 FT / = ~ / exp{-r J Bi}/ (r J B 2 ) 



\ (67) 



and F(r m , B 1: B 2 ) simplifies to 

F(r m , B i} B 2 ) = 4v^F (68) 



7.7 Ultra-relativistic case for a flat beam 

A better approximation for a flat beam at high energy can be derived by using 
Eq.(49). Since (3 2 / '(7^) < 1 and /^/^V 2 ^ > 1 it follows that B h2 > 1 but 
B\—B 2 <C 1. Eq.(49) can then be written approximately as 

F(r m ,B 1} B 2 ) = ^B 1 -B 2 f° f{A-3r m -r 3 / 2 ^l + l/r) e'^'^ 

4 /4 1 w r r m t/t^\\ dr 

For small r m and small B x — B 2 the function F simplifies to 

F(r m ,B 1 ,B 2 ) = (4 + 2 (£x - S 2 )) + ^ - B 2 ( In (4/r m ) - 1 1) (70) 

where Bi-S 2 is given by /^/(tV^ct*). For r m <10- 3 (£<3.2%) and £i-£ 2 <0.3 
{ r y>2(3 x / <j x p) the error is smaller than 8 * 10~ 3 . 

Eq.(70) shows that for decreasing r m F increases continuously, and that the gra- 
dient of the increase goes to zero with increasing energy. 
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7.8 Ultra-relativistic case for a round beam 

In this case both, Px/(l 2a xp) an d fill '{.l 2(J %)i an d therefore also B\ are small as 
compared to one (B 2 = 0) . With approximations for the Bessel functions K Q and K x 
[11] and with Eqs.(52) and (54) one obtains 

F{r m ,B 1 ,0) = v^(4 + Si(1.73 + 2 1n(Si)-8v^)) (71) 

where B 1 is given by /3X/( 7 V^). For r ™ < 10 ~ 3 ( 6 < 3 - 2% ) and B i < 0A (7 > 
SPxf&xp) the error is smaller than 3 * 1CT 3 . 

* * * 
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Appendix Al 

Since in the c. o. m. system E\ = E 2 Eq.(6) can be written as 

Ei,2 = (E 1 + E 2 )/2 

= lt {E 1 + E 2 -f3 t {f3 l E 1 cos xi + foE 2 cos X2) ) /2 

= lt (E 1 +E 2 -p 2 (E 1 +E 2 ))/2 

= (E 1 +E 2 )/(2 lt ) 

» E/ lt (Al.l) 

With Eq.(7) one gets 

J_ 1 _ c 2 (p 1 +p 2 ) 2 

7t 2 ~ (Ex + E 2 y 

E\ - c 2 pj + E 2 - c 2 p 2 2 + 2EiE 2 - 2c 2 pip 2 cos(xi + X2) 

(Si + E 2 f 

1 + 7l72 ~ /?l/?27l72COs(Xl + X2) 

(7i + 72) 2 
„ l + 7 2 -/gV(l-2x 2 ) 

~ 

2 7 2 
1 + iVx 2 
7 2 

From the Lorentz transformation 

Pii,2 = 7t(Pji,2 + (3tEi, 2 /c) 
follows with Eqs.(7) and (Al.l) 

= Pjifiht- P t (E 1 +E 2 )/(2cr rt ) 
= (2Pji,2-Pji-Pj2)/(27t) 
= ±(Pji-Pj2)/(2-f t ) 

■ Pjl-Pj2 



27*(Pji +Pi2) 

With pa = (Eq.(4)) and with Eq.(Al.2) one gets 

Pjl,2 - ± 



± 



2lt(Pjl+Pj2) 

7&>(P1 + P2) 



2lt{Pjl + Pj2) 



(A1.2) 



fp(pi+p 2 ) /, j n 2 2 2 

2(^1+^2) 

±^V / Wf/2 (Al-3) 
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Appendix A2 

We replace the variables £gi,2, a^gi 2> z pi,2, z 'pi,2i an d Ap^ 2 by the variables x@, x'p, 
zp, z'p, Ap, £, 0, and ( with help of the relations 

X/31,2 =XpT ArTf A 2/J1.2 = ^ =F 

x'^ 2 = 4 ± 0/2 =F L>; 7 C/2, 4 li2 = 4 ± C/2 =F ^7^/2, A Pl , 2 = Ap ± p 7 £/2 
and obtain 

2c(3(3 2 J 2 N 2 



Jxv(x)e- H dV* (A2.1) 



327rV>X^p 
with the condition x 2 > x„ = ^m/l 2 (Eq-(21)) and with 



H = 



4^ 2 4^ Aa% 

and 

dF* = dAsidxpdzpdApdx'pdz'i 3 d£ ! d9 d( 

The Jacobian of the transformation is 7p. In Eq.(A2.1) six of the nine integrations 
can be performed immediately and one obtains 



c^fcN 2 roo ,oo roo (rf - {(3 x D x 9/al p + (3 z bt/a%)a 2 ) 2 

327T 3 (7 2 fl (7 2 fl (7 <! (T„ i-oo i-oo i-oo 4(7? 



'xP u zP u s"p J-ooj-ooj-oo >±u h 



^(^^)-^-^K)^^C (A2.2) 



with the condition x 2 > xl- Integration over £ gives 



c(3(3 x (3 z a h N 2 roo roc 

R = 77 5/2 2 2 / / X^KX) 

lD7T t, ' 2 (7^ /3 (7 2/3 (7 s (7p J-oo J-oo 

ra 2 hf (3 x D x 9 f3 z b z (^2 (3 2 x 9 2 (3 2 ( 2 , in „ lk s 
xexp{^M^-^- + ^^) -^—-^l-\d9dC A2.3 

with the condition 9 2 + ( 2 = Ax 2 > Ax 2 m for the double integral. a h is given by Eq.(32). 
Substituting 

9 = y/pcosv, ( = yfpsinv, d9d( = dpdu/2 
one obtains (% = y / p/2) 
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c^ z a h N 2 ,00 , , rP/g^^^ 

Ota V aia ' OiaOia J) 



+ 

a z/3 ' u zf3 ' u x/3 u z/3 

cBB 3 UhN 2 r°° r 2w 
= on 5/22 2 £ / 9 / X^W expj-p^ +A 3 cos(2z/) +A 4 sin(2z/))}^p (A2.4) 

with 

4 4 2 / 1 & 2 , o\DK al(3 x (3 z D x D z 

With cos0 o = A 3 /y/A% + A\ and sin0 o = A A / ^ A\ + A| and using [11] 

7 (pA 2 ) = - e ±pMms6 d6 = — e ±pMco ^ 2d) d6 
ix Jo 2ir Jo 

where I Q is the modified Bessel function, the double integral simplifies to 

R = Qot/f^f" — / / X^{x)« x p{~P a i -pA 2 cos{2u- <f) )} dv dp 
32ir 5 / 2 a< l3 a 2 zl3 a s a p J^ m Jo 

cpp x 6 z a h N 2 roo 

= 16^7V>5^ ex P^> *P (A2.5) 



with 



64VaLV a 2 y ^ v a 2 y; 



^2 f3x(3z a h (J2.JI ^-4 n 2 n 2\ , . _ „ N 



J x(3 a z(3 u p 

With Eqs.(23) and (25) one gets (p = 4^ 2 



32^ 2 a 2 x(3 a 2 z(3 cr s a p J^/^K^ P 2J K 51 > ltX 

+ ( 3 -i-i) ln T^) ex P {-pA 1 }I ( P A 2 ) ltX dp (A2.7) 



8 2 (3 A 

With p = 4r/( / 3V), B 12 = 4A h2 /(3 2 -f 2 ), and r m = B 2 5 2 m one obtains finally Eq.(31). 
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F(T m ,Bi,B 2 ) 




Figure 2: F as a function of x m for different ratios of beam height to beam width. 
( Y = 10 4 , a xp ^ / (P x o x ) = 10" 3 , D z = 0, D;= ) 
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F(T m ,Bi3 2 ) 
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F(T m ,Bi,B 2 ) 



7 



6 







y = 500 

600 

750 

1000 






2000 

5000 

10000 

. CO 

































10~ 6 10~ 5 10~ 4 10~ 3 



Figure 4: F as a function of x m for a round beam and different energies. 

(g xP ^/(P x g x ) = io- 3 , d z = o, d;=o) 
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